SOLVING AN INVERSE OBSTACLE PROBLEM FOR 
THE WAVE EQUATION BY USING THE BOUNDARY 
CONTROL METHOD 



LAURI OKSANEN 

Abstract. We introduced in [T7] a method to locate disconti- 
nuities of a wave speed in dimension two from acoustic boundary 
measuments modelled by the hyperbolic Neumann-to-Dirichlet op- 
erator. Here we extend the method for sound hard obstacles in 
arbitrary dimension. We present numerical experiments with sim- 
ulated noisy data suggesting that the method is robust against 
measurement noise. 



1. Introduction 

Nondestructive obstacle reconstruction through wave propagation 
motivates a number of mathematical problems with several applica- 
tions such as medical and seismic imaging. There is a large body of lit- 
erature concerning obstacle detection using time harmonic waves, and 
we refer the reader to the review articles (HI [EH] and to the monograph 
|13j . Recently there has been also interest in reconstruction methods 
from acoustic measurements in the time domain P, [151 HB]- In this 
paper we present a numerical method of the latter type. We allow the 
background to be anisotropic and non-homogeneous but confine our- 
selves to the case of non-stationary acoustic waves and the scattering 
from sound-hard obstacles. 

Let M be a compact smooth manifold with smooth boundary dM 
and let g be a smooth Riemannian metric tensor on M. Let S C M mt 
be a compact set with nonempty interior and smooth boundary, and 
let \x G C°°(M) be strictly positive. We consider the following wave 
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equation on M, 

(1) d 2 u(t, x) - A g ^u(t, x) = 0, (i, x) G (0, oo) x (M \ E), 

d U)fl u(t,x) = f(t,x), (t,x) G (0, oo) x <9M, 

d u ^u(t,x) = 0, (if:, x) G (0, oo) x <9£, 

u\ t =o{x) = 0, 9 t w| t= o(x) = 0, x G M \ £, 

where A 9j/i is the weighted Laplace-Beltrami operator and is the 
normal derivative corresponding to A g tfl . That is, if we let {g^ k {x))^ k=1 
and \g(x)\ denote the inverse and determinant of g(x) in local coordi- 
nates, then we have 

A g ^u = iT 1 div(/igradw) 



E 



j,k=i 

n 



d Vytl u = /i(gradw, v) TMxT * M = ^ /j(x)u k (x)g :,k (x)^-, 

j,k=i 

where v — (z/ 1; . . . , u n ) is the exterior co-normal vector of dM normal- 
ized with respect to g, that is, J2Tk=i g^ v j v k = 1- 

Let us denote the solution of ([I]) by u^{t,x) = u{t,x). For T > 
and an open V C dM we define the operator 

Ar,r : / m- ^|(o,T)xr, / e C °°((0,T) x T). 

The Neumann-to-Dirichlet operator A<r,r models boundary measure- 
ments with acoustic sources and receivers on V. Let us assume that 
the metric tensor g and the weight function /i are known but E is 
unknown. We consider a method to locate £ from the measurements 
A Tj r. 

Let us point out that if M C W 1 , g = c{x)~ 2 dx 2 and fx(x) = c(x) n ~ 2 
where c G C°°(M) is strictly positive, then A 9i/1 = c(x) 2 A, where A is 
the Euclidean Laplacian. Thus the isotropic wave equation, 

d 2 t u - c{x) 2 Au = 0, 

is covered by the theory. The more general equation ([I]) allows for an 
anisotropic wave speed to be modelled. 

1.1. Statement of the results. Notice that the operator A 9tll with 
the domain H 2 (M) PI Hq{M) is self-adjoint on the space L 2 (M; fidVg), 
where dV g is the Riemannian volume measure of (M, g), that is, \xdV g = 
n\g\ l l 2 dx in local coordinates. We call jjdV g the measure corresponding 
to A g>fl and denote it also by V. 
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We define for a function r : dM — > M. the domain of influence with 
and without the obstacle, 

M E (r) := {x G M\ S; there is y G dM such that d^(x,y) < r(y)}, 

M(t) := {x G M; there is y G <9M such that d(x,y) < r{y)}, 

where ds is the Riemannian distance function of (M \ and d is 
that of (M,g). As (M,g) is known, we can compute the shape of the 
domain of influence M(r) for any r : dM — > M. Our main theorem is 
the following: 

Theorem 1. Let T > and /ei T C dM be open. For a function t in 

C T {T) := {r:dM^ E; r| r G C(T), < r < T, r\ dM ^ T = 0}, 

t/ie volume V(Ms(t)) can be computed from A2T,r by solving a sequence 
of linear equations on L 2 ((0,T) x T). Moreover, 

(2) M(r) mt n S m * ^0 and on/y if V(M E (r)) < V(M(r)). 

Theorem [T] allows us to probe the obstacle with the known domains 
of influence M(r), r G CV(T). We will illustrate this probing method 
in Section [3] via numerical experiments in the two dimensional case. 

In Section [2] we give a proof of Theorem [T] that is based on ideas from 
the boundary control method. By using the boundary control method, 
a smooth wave speed can be fully reconstructed from the Neumann- 
to-Dirichlet operator. This uniqueness result is by Belishev p] in the 
isotropic case and by Belishev and Kurylev [3] in the anisotropic case. 
We refer to the monograph [12] and to the review article [2] for fur- 
ther details on the boundary control method. The boundary control 
method depends on Tataru's hyperbolic unique continuation result |2Uj . 
whence it is expected to have only logarithmic type stability. Also our 
result depends on [20], however, we overcome the ill-posedness of the 
reconstruction problem by regularizing it carefully. The regularization 
stategy is a modification of that in [1], and the iterative time-reversal 
control method introduced there can be adapted to give an efficient 
implementation of our method. 

2. Proof of the main theorem 

We begin by showing that the volumes V(Ms(r)), r G CV(T), can 
be computed from A 2 T,r by solving a sequence of linear equations on 
L 2 ((0, T) x dM). Our proof relies on general results from regularization 
theory and it can also be adapted to simplify the arguments in [T£] . 
We define the operator 

K '.= JA.2T,r®2T — RA r pjpRJQ2T > 
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where G2T is the extension by zero from (0, T) to (0, 2T), R is the time 
reversal on (0,T), that is Rf(t) := f (T — t) , and 

1 r 2T ^ 

Jf(t) := - J f(s)ds, f G L 2 (0,2T), t G (0,T). 

We recall if is a compact operator on L 2 ((0,T) x T) since, see [2~T] , 

A r ,r : L 2 ((0,T) xF)-> F 2/3 ((0,T) x T). 

Let / G C£°((0,T) x T) and let G C°°(M \ S). Moreover, let 
t G (0, 00) and integrate by parts 

(3) ^(u f (t),(f)) L 2^ M \ s . dV) = (A g)IJ ,u f (t),(f)) L 2 {M \z. dV) 

= -(gradn / (t),grad (j)) L ^ M \r,;dV) + {d u ^u f {t),<j)) L 2 {dM . dSg)} 

where dS g denotes the Riemannian surface measure on (<9M, g). Notice 
that the boundary term on <9£ vanish as u* satisfies the homogeneous 
Neumann boundary condition there. 

In particular, for f,he C°°((0,T) x T), t G (0,T) and s G (0,2T), 

{d 2 t -dl){uf{t),u\s)) L , {M ^, dV) 

= (f(t),^2T,rh(s)) L 2 (dM . dSg) - (Ar,r/(t),/i(s))L2 (aM;d5g) . 

By solving this wave equation with vanishing initial conditions at t = 
and noticing that A^ r = RA T r R, we get the Blagovescenskii's identity 

(4) {w{T),u h {T)) L i(M\Y,-dV) = (/, Kh) L 2(( 0:T ) xr . dtmSg ), 

that holds for all /, /1 G L 2 ((0,T) x T) by continuity of if and density 
of smooth functions in L 2 . The identity Q originates from [5]. 
Moreover, by letting = 1 identically in ([3]), we get 

(5) d 2 (u f (t), l) L 2(M\X;dV) = (d^U f (t), l) L 2(dM;dS g )- 

Notice that this identity does not hold if satisfies the homogeneous 
Dirichlet boundary condition on <9£, instead of the Neumann one. This 
is why our method does not extend to detection of sound soft obstacles 
in a straightforward manner. We get the indentity 

(6) ( uf (T), l)L 2 (M\T,;dV) = (f,b) L 2(( 0tT ) xr . dmdSg ), 

where b(t, x) = T — t, by solving the ordinary differential equation ^ 
with vanishing initial conditions at t — 0. 
Let r G CV(r) and let us define the set 

S T := {(t,x) G [0,7] xf; t G [T — r(x),7]}. 

We define the operator 

W T f:=u f (7), W T : L 2 (S T ) -> L 2 (M\£). 
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It follows from [T4] that W T is compact. Moreover, we may consider a 
restriction of K, 

K T f = Kf\ Sr , K T :L 2 (S r )^L\S T ). 

Then the equations Q and (JsJ) yield that on L 2 (S T ) 

(7) W*W T = K T , W*l = b. 
Let us now consider the control equation, 

(8) W T f = l, for fEL 2 (S T ). 

We have supp(W / T /) C M s (r) since the wave equation Q has finite 
speed of propagation. Moreover, it can be shown using Tataru's unique 
continuation [20] that the inclusion 

(9) {W T f; f G L 2 (S T )} C L 2 (M s (r)), 

is dense, see the appendix below. In particular, if there is a least squares 
solution f to ([8| then W T fo = 1m e (t)- However, as W T is compact, the 
range of W T is a proper dense subset of L 2 (M 2 (r)) and ^ may fail to 
have a least squares solution. Nonetheless, it is instructive to consider 
first the case where ^ has a least squares solution. Then the least 
squares solution of minimal norm /o is given by the pseudo inverse, see 
e.g. Th. 2.6], 

f = w\\ = (w;w T yw;i = k\i, 

and we can compute the volume V r (Ms(r)) from the boundary data 
A2t,a by the formula 

V(Mv(t)) = (1m e (t),1)l2(M\IW) = (W T W}1, l)L2(M\S;dy) 

= (Klb, b) L 2( ST . dmdSg ). 

The standard technique to remedy the nonexistence of a least squares 
solution to a linear equation is to use a regularization method. As W T 
is compact and we have the information Q at our disposal, there are 
several ways to regularize that are available to us. For example, we 
could use a regularization by projection [HI Section 3.3] or a regular- 
ization based on a spectral approximation of the inverse [91 Th. 4.1]. 
Here we will consider only the classical Tikhonov regularization, 

(10) f a := (W*W T + a^W*! = (K T + a)~\ a > 0. 

We have the following abstract lemma. 

Lemma 1. Suppose that X and Y are Hilbert spaces. Let y G Y and 
let A : X — >■ Y be a bounded linear operator with the range R(A). 
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Then Ax a — > Py as a —> 0, where x a = (A* A + a)~ 1 A*y, a > 0, and 
P : Y — > R(A) is the orthogonal projection. 

Proof. Notice that for all x G X 

\\Ax - yf = \\Ax - Pyf + ||(1 - P)y\\ 2 . 
By [HI Th. 5.1] we know that x a is the unique minimizer of 

|| Ax — y\\ 2 + a \\x\\ . 
Let e > and let x e G X satisfy \\Ax e — Py\\ 2 < e. Then 
\\Ax a - Pyf = \\Ax a - yf - ||(1 - P)y\\ 2 

< \\Ax a -yf + a\\x a \\ - ||(l-P)y|| 2 

< \\Ax e - yf + a Wx'W - ||(1 - P)yf 

= \\Ax e - Pyf + a \\x e \\ <e + a \\x e \\ < 2e, 
for a < e/ \\x e \\. □ 



By the density Q we have that R(W T ) = L 2 (M s (r)). Thus the 
above lemma implies that Wf a — > 1m e (t) i n L 2 (M \ E) as a tends to 
zero. In particular, we may compute the volume V(Ms(t)) from the 
boundary data A 2 t,a by the formula 

(11) V(Mx(t))= \im((K T + a)-%b) LHST . 4mdSg y 

a— >0+ 

Lemma 2. Let T > 0, T C dM be open and let r G C T (T). Then 
M(r) int n E int ^0 «/ and on/y i/ K(M s (r)) < V(M(r)). 

Proof. Notice that ds(x,y) > d(x,y) for any x,y G M \ E. Hence 
Me(t) C M(r). Morever, Mj(r) D E = by definition. In particular, 
if the open set M(r) int fl E mt is nonempty, then 

V(Me(t)) < V{M{t)\E) 

< V(M{t) \ E) + V(M(T) int n E int ) < V{M{t)). 

Thus we have shown the implication from left to right in ([2]). 

Let us now suppose that V(M E (r)) < V(M(t)). Then M(t)\M^(t) 
is not a null set (that is, a set of measure 0). But dM(r) is a null set 
|18j . whence there is x G M(r) mt \ Me(t). Thus there is y G and 
a path 7 : [0,£] — > M from y to x such that the length of 7 satisfies 
Kt) — r (l/)- The path 7 intersects E since otherwise we would have 
x G M s (r). Let 2 G E n 7([0,£]). Then z G M(r) int since x G M(r) int , 
and there is a neighborhood C/ C M(r) mt of z such that U fl E int 7^ 0. 
Hence also M(r) int n E int ^0. □ 
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Theorem [I] follows from the formula (11) and Lemma [21 



3. Numerical results 

3.1. Simulation of the data. In all our numerical examples (M,g) 
is the two-dimensional unit square with the Euclidean metric, that is, 

M = [0,1] 2 , g = (dx 1 ) 2 + (dx 2 ) 2 . 

Moreover, T = 1 and the accessible part of the boundary T is the 
bottom edge of M, 

T = {(^,0) eM;x l G (0,1)}. 

For computation of the Dirichlet-to-Neumann map we discretize in 
space by using finite elements, and solve the resulting system of ordi- 
nary differential equations by a backward differentiation formula (BDF). 
To be very specific, we use the commercial Comsol solver with quadratic 
Lagrange elements and BDF time-stepping with maximum order of 2. 
Both the maximum element size and time step size are set to the con- 
stant value h = 0.0025. 

We discretize the measurement Ar^rf, f £ L 2 ((0, T) x T), by taking 
the point values on the uniform grid of temporal points tj e [0, 2T], 
j = 1,2, ...,N t , and spatial points Xk € T, k = 1,2, . . . , N x , where 
N x = 20 and N t = 800. The higher precision in time reflects the fact 
that a measurement of this type can realized by using N x receivers (e.g. 
microphones) with the sampling rate h. 

We model noisy measurements by adding white Gaussian noise to 
the signal 

\f(j,k):=A r , T f(t j ,x k ), j = 1, 2, . . . , N t , k = 1,2, . . . , N x . 

To be very specific, we use the Matlab function awgn both to measure 
the power of the signal Xf and to add noise with specified signal-to- 
noise ratio (SNR). We have used signal-to-noise ratios 14dB and 7dB 
corresponding to 4% and 20% noise power levels. 

3.2. Solving the control equation. The operator K T is self-adjoint 
and positive-semidefinite by ([7]), whence K T + a positive-definite for 
a > 0. We solve the Tikhonov regularized control equation 

(12) (K T + a)f = b 

by using the conjugate gradient (CG) method on a finite dimensional 
subspace C T C L 2 (S T ) that we will define below. We have used the 
initial value / = in all our CG iterations. 
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Figure 1. Convergence of the reconstructed volume V(M$(T r )) 
as a function of N cg , i.e. the number of the conjugate gradient 
steps. Noiseless case. Left: r = 1/10; right: r = 1/2. 



We denote by C T the Voronoi cell corresponding to the measure- 
ment point Xk-, k = 1, 2, . . . , N x , that is, 

T k := {x G T; \x - x k \ < \x - x t \, I = 1, 2, . . . , A^}. 

Moreover, we denote by C the space of piecewise constant sources / 
that can be represented as a linear combination of the functions 

(13) f k (t, x) := l m {t)l Tk (x), k — 1,2,..., N x , 

and their time translations by an integer multiple of h. Finally, we 
define 

C T := {/ G C- supp(/) C S T }, r G C T (T). 

As the wave equation ([T]) is invariant with respect to translations in 
time, we can compute A/ for arbitrary / G C T and r G CV(r) if we are 
given the measurements 

X fk , k = l,2,...,N x . 

To summarize, we employ N x = 20 measurements that can be realized 
by using receivers with the sampling rate h = 0.0025. 



3.3. Regularization and calibration. As the control equation (12) 
may be ill-posed for a = 0, we terminate the CG iteration early after 
N cg steps. This amounts to regularization of the problem [10J. To 
calibrate the method we probed the empty space case, £ = 0, with 
half-spaces. That is, we chose the profile function r G CV(r) to be of 
the form, 

T r (x):=r, xeT, r G [1/10, 1/2]. 

In this case, the CG iteration essentially converges after 10 steps even 
when not using the Tikhonov regularization, that is, a = 0, see Fig- 
ure [Tj For this reason, we have chosen N cg = 10 in all our further 
simulations. 
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Figure 2. Reconstructed volumes V(M${r r )) as a function of 
r compared to the real volume (dotted red). The two recon- 
structions (solid blue and dashed blue) correspond to two differ- 
ent realizations of noise. From left, 1st: noiseless, a = 0; 2nd: 
SNR = UdB, a = 0; 3rd: SNR = 7dB, a = 0; 4th: SNR = 7dB, 
a = 1(T 3 . 






0.1 0.3 0.5 



Figure 3. Reconstructed volume differences (14) with r = r r as 



a function of r compared to the real difference (dotted red) in the 
case of the disk shaped obstacle £ = E . The two reconstructions 
(solid blue and dashed blue) correspond to two different realiza- 
tions of noise. From left, 1st: noiseless, a = 0; 2nd: SNR = lAdB, 
a = 0; 3rd: SNR = 7dB, a = 10" 3 . 



In addition to the empty space case, we have experimented with the 
disk and the square shaped obstacles defined as follows: S D is the disk 
with radius 3/10 and center p := (1/2, 1/2) and E is the square with 
side length 0.424, center p and axes rotated by 7r/4 with respect to the 
axes of M, see Figure |4| 

It is not clear to us, why the method underestimates the volume 
V(M^(r r )), see Figure [2] (leftmost plot). One possibility is that we 
using too few spatial basis functions, however, the smallness of N x is 
motivated by applications. Moreover, the underestimation is system- 
atic and is canceled when considering the volume differences, 



(14) 



V(M s (t))-V(M (t)), 
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Figure 4. Top row from left, 1st: The disk shaped obstacle 
X = S . 2nd: The largest region on which the absence of the 
obstacle can be concluded when probing with disk shaped do- 
mains of influence H-£ o (T). 3rd: A reconstruction of H-£, o (T) in 
the noiseless case. The threshold e = 5/10 4 . Bottom row: The 
square shaped obstacle S = S , H^ (T) and a reconstruction of 
Hx (T). The same parameter values are used for both the recon- 
structions. 



see Figure [3j In terms of applications, this means that we should 
calibrate the method in a known background before probing a region 
that possibly contains an obstacle. 

According to our experiments the method reconstructs volumes re- 
liably when SNR = 14dB and we regularize only through the early 
termination of the CG iteration. When SNR = 7dB and a = 0, a re- 
construction can be seriously disrupted even in the empty space case. 
After introducing Tikhonov regularization with a = 10~ 3 , the effect of 
noise vanishes but a large systematic error appears, see Figure [2] (the 
two rightmost plots). We see that considering the volume differences 
(14) becomes even more essential when a > 0. 



3.4. Probing with disk shaped domains of influence. We will 

now describe our experiments concerning reconstruction of the shape 
of an obstacle. To this purpose, we chose the profile function r G Ct{T) 
to be of the form, 



rH(x) 



r-\x-y\, x G T, y G T, r G [1/10,1/2]. 



Then M(r^) = B(y,r) H M, that is, the intersection of M and the 
closed disk of radius r centered at y. Probing with disks has been 
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Figure 5. Comparison of a reconstruction Hs(T) with the the- 
oretical best possible reconstruction H^(T). Erroneous pixels in 
H S (T) \ Hx(T) are drawn in white and H S (T) \ H^(T) in black. 
Gray pixels are reconstructed correctly. Top row: Noiseless mea- 
surements with the threshold e = 5/10 4 , left: E = E ; right: 
E = E . Bottom row: SNR = UdB, E = E D and e = 4/10 3 . 
The two reconstructions correspond to two different realizations 
of noise. 

considered in the context of electrical impedance tomography in [TT] 
and our numerical results are comparable to the results therein. 

Analogously to [11] and [17], let us define the largest region H-ziT) 
on which we can conclude the absence of obstacles by probing with the 
sets B(y, r) D M, y G T, r G (0, T]. We denote 

R T (y) : = sup{r G (0,T] ; B(y,r) n E int = 0} 

= sup{r G (0,T] ; V{M s (t?)) = V{M{ T y))}, 

and define 

H E {r):=\J{B(y,R T {y))nM). 

Let us describe next how we approximate Rt{u) when computing 
with finite precision. Let e > 0, N r G N and let r/ G [0, T], I = 
1,2,..., N r , be a uniform grid of points. We denote 

L(e, N r ) := max{Z = 1, 2, . . . , N r ; 7(Af s (r»)) - V(M 9 (t*)) > -e}, 

and define the approximation r T (y; e, N r ) = r L ^ e>Nr ) of Rr(y)- We have 
used the threshold e = 5/10 4 in noiseless cases and e = 4/10 3 when 
SNR = lAdB. According to our numerical experiments the method 
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reconstructs H^(T) reliably when using these values of e and N r = 500, 
see Figure |5j where a white pixel means that the center point of the 
pixel is erroneously identified to be in H^(T) (false positive) and a 
black pixel means erroneously identification of not being in H^(T) (false 
negative) . 

Computationally the shape reconstruction amounts to solving a large 
number of independent systems of linear equations by running a few 
number of CG steps for each of them. Our implementation with pa- 
rameters as above and r(s restricted in [1/10, 1/2] led to 4020 systems 
with the number of unknowns varying between 30 and 1000. The run 
time for the full reconstruction on a single processor was about 10 min- 
utes, however, as the systems are independent, the method allows for 
an efficient parallel implementation. 

Appendix: Approximate controllability 



In this section we show that the inclusion (|9j) is dense, that is we 
prove the following lemma. 

Lemma 3. Let T > 0, let F C dM be open and let r G CV(r). Then 

(15) {u f (T); f G C™(S T )} 

is dense in L 2 (Me(t)). 

A density result of this type is called approximate controllability in 
the control theoretic literature. To our knowledge, Lemma [3] is proved 
previously only in the case of a constant function r, see e.g. [121 Th. 
3.10]. We will give a proof in the general case r G Ct(T) by reducing 
it to the constant function case. To simplify the notation we consider 
only the case E = 0, since the general case follows by replacing M by 
M \ E in the proofs below. 

Lemma 4. Let T > 0, J G N, let Tj C dM be open and let hj G C T (Tj) 
for j = 1, 2, . . . , J. We define 



(16) h J {y) :- 



max{hj(y);j satisfies Tj 3 y}, ye \J j=1 Tj, 
0, otherwise. 



If for all j = I,..., J the functions (15) for r = hj are dense in 
L 2 (M(hj)) , then the functions (15) for r = h J are dense in L 2 (M(h J )). 



Proof. Notice that dM C M(r) if r(y) > for all y G dM. Abusing the 
notation slightly, we will consider M{r) as a subset of M mt . This does 
not affect the density since dM is a null set. We denote T J := [L = i 
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and have 

M(h J ) = {x e M int ; there is y G f J s.t. d(x,y) < h{y)} 
j 

= (J{x G M int ; there is y G Tj s.t. d(x,y) < hj(y)} 

3=1 
j 

= U M ( r i^i)- 

3=1 

We will now prove the density by induction with respect to J. The 
case J = 1 is trivial. Let us denote M := M{h J ) and M\ := M{h J+ i). 
Let ip G L 2 (M U Mi). By induction hypothesis there is a sequence of 
smooth functions (fk)^=i supported in S h j such that 

u f °(T) -> l Mo V- 

Moreover, there is a sequence of smooth functions (fDkLi supported 
in Sh J+1 such that 

u f *{T) ^ImM~ W)- 

Thus 

U fk+fk(T) -> l Mo (l - + ImiV' = (1m \Mi + ImJ^ 

= Imoum^ = ip. 

Moreover, + fl is supported in S^j U Sh J+1 C . □ 



Proof of Lemma^ Let ^ G L 2 (M(r)) and e > 0. There is a simple 
function 

j 

Mj/) = Z} T iMi/). 

3=1 

where J G N, 3} G (0,T) and Lj C V are open and disjoint, such that 
t < h e + e almost everywhere on F and h e < r on T, see e.g. [18| 
proof of Lem. 4.2]. We denote hj := I}lf. and define r e = /i J as the 
maximum (16). By the construction r e < r and r e > h e . 

The functions (15) for r = hj are dense in L 2 (M(hj)) by [121 proof 
of Th. 3.10]. Lemma [4] implies that there is a smooth function / 
supported in S Te C S T such that 

||Uf(r E )^ - m/ ( T )||l2 (m) < e- 
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Thus 
(17) 



\\^-u f (T)\ 



L 2 (M) 



< e + 



i/j 2 dV. 



'M(r)\M(r.) 

We have V(M(r e )) V(M{t)) as e 0, see p2]. Thus the second 
term in (17) tends to zero as e — > 0. □ 
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